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Abstract 

We report on a high statistics numerical study of the crystalline random 
surface model with extrinsic curvature on lattices of up to 64 2 points. The 
critical exponents at the crumpling transition are determined by a number 
of methods all of which are shown to agree within estimated errors. The 
correlation length exponent is found to be v = 0.71(5) from the tangent- 
tangent correlation function whereas we find v = 0.73(6) by assuming finite 
size scaling of the specific heat peak and hyperscaling. These results imply a 
specific heat exponent a = 0.58(10); this is a good fit to the specific heat on a 
64 2 lattice with a x 2 P er degree of freedom of 1.7 although the best direct fit to 
the specific heat data yields a much lower value of a. Our measurements of the 
normal-normal correlation functions suggest that the model in the crumpled 
phase is described by an effective field theory which deviates from a free field 
theory only by super-renormalizable interactions. 



e-mail: jfw@thphys.ox.ac.uk 



1 Introduction 



The crystalline random surface model has attracted considerable interest over the 
past few years. Although such models have obvious applications in condensed matter 
physics IJ their role in the physics of two-dimensional quantum gravity is perhaps 
the most intriguing. For some time the nature of two-dimensional quantum gravity 
coupled to conformal matter with central charge c < 1 has been reasonably well 
understood || however this understanding fails when c > 1. Roughly speaking 
when a small number of matter species are present the physics is a small perturbation 
from that when there is no matter; when the number of matter species exceeds a 
certain value the self coupling of gravity through intermediary matter fields becomes 
strong and the method of calculation fails. (Whether the physics changes abruptly 
is not yet properly understood.) 

It is quite remarkable that there is an (apparently) exact lattice formulation 
of two-dimensional quantum gravity problems. In this formulation the functional 
integral over the metric is replaced by summing over all triangulations of the lattice 
H; an d matter fields are provided by placing a spin system on the lattice. 
At its critical point the model is the same as conformal matter (with c equal to 
that of the spin system at its fixed lattice critical point) interacting with quantum 
gravity; this is known by direct calculation in the case of the Ising model for example 
JF|. In this approach the simplest way to obtain integer c = D is to put D free 
(that is, not self-interacting - they do of course interact through their coupling to 
the dynamical triangulation), massless bosonic fields on the lattice. The model is 
then mathematically equivalent to one of a random surface, which is allowed to 
self- intersect, embedded in D dimensions. Unfortunately, it turns out that such 
models do not lead to a smooth continuum limit because there are no correlations 
between normals of the surface and the continuum string tension is always infinite 
[|J. To overcome this problem it was proposed to add an extrinsic curvature term 
to the action ||; this introduces a correlation between normals and so might lead 
to a smoother surface with a finite string tension in the continuum limit. The 
resulting model on a fixed triangulation, usually called the crystalline surface model, 
has been shown to have a second order phase transition, as have similar models 
Tl~i 12| , |i~3] , 14| . Of course, since this is no longer a free field theory, it is clear 



neither that the effective field theory governing the critical point is a conformal field 
theory nor, if it is, whether it actually has c > 1. If this critical point is conformal 
with c > 1 then the corresponding dynamically triangulated model might lead to a 
consistent theory of quantum gravity interacting with c > 1 conformal matter. 

There are several recent papers which review the status of random surface models 
per se 0, [H| and various approaches to two-dimensional quantum gravity [|1^, |I7 



The purpose of this paper is to report on a high statistics numerical study of one 
model in D = 3 dimensions in the region of its critical point; preliminary results 
which did not use the full data set analyzed here were presented in IjTBl . This paper 
is organized as follows. In sections 2 and 3 we describe the model and our simulation 
methods, sections 4, 5 and 6 deal with the results and in section 7 we discuss our 
conclusions. 
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2 The model 



The crystalline surface model approximates a continuum surface by a (regular) tri- 
angulation with N vertices labelled i = 1, . . . JV; each vertex has a location on the 
lattice denoted by its two-component intrinsic coordinates £(z) and a location in 
the D-dimensional euclidean embedding space X(z). We will sometimes refer to the 
components of X as where fi = 1, . . . D. In numerical work it is much easier to 
use closed surfaces and our choice of lattice boundary conditions (periodic in y and 
helical in x |12[) make the surface a torus. 



There are several ways of assigning an action to this system JT0|, [11], |18|] and 
for some time it appeared that they did not all have the same critical properties. 
However, recent work |14] has demonstrated that the analysis of simulations is com- 



plicated by very long autocorrelation times and that in fact all the models have 
similar critical behaviour. In this paper we work with the model whose action is 
given by 

S = (3 2 S G + KS ec . (1) 
where (3 and k are coupling constants. Sq is gaussian and given by 

Sg = \ £(X(z)-X(j)) 2 (2) 

Z <ij> 

where < ij > denotes the link joining neighbouring points labelled by i and j. S ec 
denotes the extrinsic curvature action and is given by 

See = ~ £ A a ■ n A ' (3) 

<ij> 

where A, A' are the triangles sharing the link < ij >, n A is the normal vector of 
A, and n A the correspond unit vector [10|. The partition function is given by 

Z = Ufd D X( i )5 D (j:X( i ))e- s , (4) 

where the delta function has been included to remove the translational zero mode. 
Expectation values are given by 

1 N r / N \ 

(O) = % II J 8 D (E X(z) J 0({X}) e- s . (5) 

By exploiting the invariance of S ec under the global rescaling X — > AX it is easy to 
derive the results 

(N-1)D 

( S g) = — 2/35 — (6) 
(5S> - (So) 2 = (4)c = (7) 
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and so on; note that (3 can be scaled away entirely but we will retain it for later 
convenience. 

At k = this model is exactly soluble and describes a surface whose mean square 
radius 

1 N 

(x 2 ) = -yExW-xw (8) 

i=l 

behaves as 

(X 2 )~logiV. (9) 
Defining a "Hausdorff" dimension du through 

(X 2 ) ~ N^h (10) 

the behaviour at k = corresponds to du = oo and the surface is so highly crumpled 
that it is space filling in arbitrary dimension D; this behaviour persists at small k 
and this phase is usually called the "crumpled phase" . The extrinsic curvature term, 
*S ec , acts to suppress sharp bends in the surface and so as k increases we expect the 
surface to become smoother. It is now well-established JTO, [TTJ, [T2], |HJ that there is 
a second order phase transition at some critical value, k c , of the coupling. There 



are reasons for believing that at k = k c the value of du drops to 4 or less [11], [12 1; a s 
k — ► oo the surface becomes smoother and we expect du — > 2. In fact the numerical 
evidence [^2| indicates that dn goes straight to 2 at the critical point. 



Although considerable effort has gone into determining the nature of the critical 
point and measuring the critical exponents there are still some apparent discrep- 
ancies. The purpose of this work was to gather very high statistics on lattices of 
moderate size in order to see if the discrepancies could be resolved. As a byproduct, 
we have also been able to measure some quantities for the first time. 



3 The simulation 

There are some advantages to simulating these models using the Langevin method 
with Fourier acceleration. This is largely because the action depends only upon 
derivatives of the Xs and not on the fields themselves; it has been observed |12j that 
a Fourier accelerated simulation gives much shorter autocorrelation times. Although 
this is partly offset by the extra computer time needed for the Fourier transforms 
there is still a substantial gain. However the method is not perfect and critical 
slowing down is still observed near the phase transition. The method of Fourier 



acceleration used in this work is identical to that described in |~2] but another 
aspect of the algorithm has been changed. 

In a Langevin simulation, new configurations are generated by discretizing the 
Langevin time r in steps of 5t; at an update the change in X is given by 

SX,(i) = -*TQx-ffi + V(h r)Vfr (11) 
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where 77 is a gaussian distributed random number satisfying (r](i, t)t)(J, r')) = 2Sij8 TT i 
The derivative of S ec depends upon 



(Ma 1 / dn. 




&A Ua-^ (12) 



dX M (z) |n A | V^(0 V dx ^ 

which is non-zero if % is one of the vertices of A and zero otherwise. If A has 
a very small area |iia| then it is possible that a very large update for X^ii) will 
be generated; there is no reason why the expression in brackets in fll2|) should be 
particularly small when |iia| is small. Thus configurations of the surface which have 
one or more triangles of very small area are prone to generate large spikes at the next 
update. This problem exists for any finite St and any non-zero k and is an artefact 
of discretizing the Langevin equation. To deal with the difficulty it is necessary to 
modify the simulation procedure so as to suppress the contribution of triangles of 
very small area. This may be done by introducing a small area cut-off; the resulting 
procedure is rather like a hybrid algorithm. However this is computationally very 
inefficient to implement and we have developed a different method. 
We work with the modified extrinsic curvature action 



s ec = ~ J2 n A ■ nA'exp \-oSt ( 7— t + 7-7-7 J J 

<if> V vl n Al |n A |yy 



(13) 



<ij> 

where a is some constant so that the total action for the simulation is given by 

S S ym = (3 2 S G + nS ec . (14) 

When St —* this reduces (at least naively) to the original action. On the other 
hand at finite St the contribution of any triangle of area much less than oSt to the 
Langevin update is heavily suppressed. The exponential suppression wins over any 
power law from fll2|) and provides a cut-off which is almost a theta function but 
has the virtue of being differentiable; thus S ec can be simulated by a completely 
standard Langevin method. (Note that this procedure does not suppress the cre- 
ation of small triangles, it merely prevents them generating spikes.) Working with 
S ec is also very efficient computationally; compared to the original action S ec very 
little extra calculation is needed. Simple analysis shows that a is dimensionless; 
in principal we could choose it to take any value but its presence can be exploited 
further. A Langevin simulation with finite St always renormalizes the couplings of 
the bare action and in general introduces other operators as well ||20|| . This has the 
consequence that the X configurations that we generate do not actually yield the 
expectation value (|6]) correctly because f3 is renormalized; this trivial renormaliza- 
tion can be computed by comparing the measured and theoretical values of (So)- 
However, even after taking account of the renormalization of (3, higher moments 
such as (0) do not come out as predicted because other operators are induced by the 
discretization. We use the freedom to tune a to arrange that the measured ratio 

P= {J ^^(S 2 G ) C /(S G ) 2 (15) 
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is close to the predicted value of 1; that is to say we can remove at least some of 
the effects induced by the discretization of the Langevin time. We found that it is 
easy to tune a so that the residual discrepancy in p is of order 1%. The value of a 
needed is approximately independent of St, as it should be since [a] = 0, at least 
over the range 5r = 0.001 to 0.004 which brackets the value of 5t = 0.002 used 
in our simulations as is shown in fig. 1 . That the procedure is consistent is shown 
in fig. 2 ; the value of a deduced by tuning on a 16 2 lattice works just as well on 
a 64 2 lattice. Higher moments of Sq are still affected by discretization errors; for 
example the connected third moment comes out about 50% too large (although it 
should be realized that it is very difficult to measure this double-subtracted quantity 
accurately) . 

Despite the improvements due to the Fourier acceleration it is still necessary to 
acquire very high statistics for reliable measurements. For example, in the case of 
the specific heat (defined in section 5) the total data set consisting of T sweeps is 
split into K bins; we then compute the specific heats, C from the total data set, 
and Cfc (k — 1, . . . K) from each bin separately. We define K bin to be the largest 
value of K for which 



1 K 



< 0.01 (16) 



and then 

T bin = TqK^I. (17) 

Subsequent analysis is then done using bins of size T bin ; in this way we ensure that 
the statistical errors in the measured specific heat are larger than the thermalization 
errors. The values of T bin found, measured in numbers of second order Runge-Kutta 
iterations of our algorithm, are shown in fig. 3; as can be seen Tu n appears to scale 
roughly like y/N which is consistent with the behaviour found in |jl2| . We have 
run our simulations for at least lOT^ iterations except in the case of crucial k 
values where we have run for much longer; in every case the first 2Tu n iterations 
are discarded. Thus we have accumulated about 10 8 iterations of the algorithm in 
total. 



4 Tangent- Tangent Correlation Functions 

Much of the information about the long distance properties of the system is con- 
tained in the tangent-tangent correlation functions. Let 

^(0=^(6,6) = x(a + U2)-x(£i,6) (is) 

t 2 (l)=t 2 (£i,&) = X(£i + 1,6 + 1)-X(£i,6) (19) 
then the correlation functions Gij are defined by 

^(0 = ^E( t .(^6)-t,(o,6)) (20) 

fa 
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On a periodic lattice it follows that 



E Gn(0 = 4E<(E 6) I " ti(0, 6)> = (21) 



N f , < 



and therefore that Gn^), which is clearly positive for small enough £, must become 
negative at large £. Thus there must be some intermediate value, £o, for which 
Cn(£o) = 0. As shown in JTT], |T2| the finite size scaling behaviour of £o with 
L = \/~N is related to the Hausdorf dimension of the surface: 



1. £ — > constant as L — > oo; this is the crumpled phase and has djj = oo. If 
there is only one mass scale then £ oc mT x where m is the mass gap. 

2. £o — > L 1 ^ 11 , < n < 1, as L —>■ oo; then dn = 4(2 — n)^ 1 < 4 and the system 
is not crumpled. 

In this work we are examining the crumpled phase so 1) is the relevant scenario. 
When L is finite, £o wm saturate as k is increased towards the critical point because 
it cannot in any case exceed L/4 and so we expect to see finite size effects. The 
values of £ extracted from G n on 16 2 , 32 2 and 64 2 lattices are plotted in fig. 4 and 
show these effects clearly. In order to extract the true value of m near the critical 



point it is necessary to have a model for the finite size effects. Previous work [T2|, |H| 
has shown that the data for Gn at long distances in the crumpled phase is quite 
well accounted for by supposing that it is described by the effective Lagrangian 

tfree = ^ (-^( K ) 2 X • V 2 X + X ■ (vfx) (22) 

which leads to the momentum space correlation function 

Gij {p) = Up) ■ t,(-p)> = M-) p2(p /^ m{K)2y (23) 

Using this asumption we can remove the finite size effects by fitting (^) to the data 
in the region where Gii(£) < 0. But, how good is the assumption? 

In fig. 5 we show a sequence of measured Gn at different k values on a 64 2 
lattice together with the fit described above. As can be seen, at large distances 
the data is well fitted by (|23|) right up to the phase transition. The last picture 
at k — 0.82 shows completely different behaviour; the surface is no longer in the 
crumpled phase. This phenomenon provides us with an upper bound on k c which 
can be used to constrain the fits for critical exponents. No error bars are shown on 
the data. The errors are smaller than the symbols on the figure but the correlation 
function at successive £ values is highly correlated; the fits of ( |2"3"D were made by 
minimizing the x 2 function obtained using a correlated error analysis and SVD. 
However we do not actually know the true probability distribution of this data and 
so the errors on the values of A and m were obtained using a bootstrap procedure; 
all the values so obtained are listed in Table 1 together with the actual number of 
degrees of freedom (as determined by the SVD) contributing to the fits. In fact the 
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Table 1: C e /f parameters deduced from correlation functions. 



K 


m 


A 


d.o.f. 


m 


9 


0.7400 


.348(2) 


8.95(7) 


15 


- 


- 


0.7500 


.313(2) 


8.64(10) 


16 


— 


- 


0.7600 


.273(2) 


8.02(8) 


17 


- 


5.8 


0.7700 


.241(2) 


7.66(5) 


18 


- 


4.7 


0.7750 


.220(2) 


7.42(6) 


17 


- 


- 


0.7775 


.203(2) 


7.21(8) 


18 


.204(4) 


- 


0.7800 


.201(2) 


7.28(6) 


18 


.203(3) 


3.7 


0.7825 


.192(3) 


7.10(8) 


19 


.192(3) 


- 


0.7850 


.174(2) 


6.84(5) 


19 


.174(2) 


- 


0.7875 


.167(2) 


6.76(6) 


19 


.165(3) 


- 


0.7900 


.155(3) 


6.59(7) 


19 


.153(4) 


3.0 


0.7925 


.144(2) 


6.47(6) 


19 


.143(3) 




0.7950 


.130(2) 


6.29(3) 


19 


.129(2) 




0.8000 


.110(3) 


6.06(8) 


19 


.108(5) 


1.8 


0.8050 


.078(5) 


5.70(10) 


18 


.076(7) 




0.8075 


.057(4) 


5.47(7) 


19 


.055(6) 




0.8100 


.039(7) 


5.36(8) 


18 


.037(12) 


1.1 


0.8110 


.030(4) 


5.33(6) 


19 


.029(11) 




0.8125 


.011(11) 


5.27(3) 


19 


.012(10) 





data is of high quality and, as can be seen from Table 1, the number of contributing 
degrees of freedom is large so our procedure should be reliable. There is however 
one possible source of systematic error. At short enough distances the two point 
function is not well described by ( p3| ) and so it is necessary to decide the minimum 
£ at which to attempt the fit. Obviously there are different ways of doing this and 
it is not clear that they will give the same results. To check on this we tried some 
other procedures than fitting to the region where Gn(£) < 0. The fh column in 
Table 1 shows the results of fitting in the range £ > 10 for all k values. As can be 
seen the fh values are always very close to the m values but the estimated statistical 
errors are somewhat larger (probably because the number of contributing degrees of 
freedom is generally smaller). 

The measured values of m for 64 2 lattices are plotted in fig. 6 . We have much 



higher statistics than in |TB[ and, crucially, reliable measurements for k above 0.8. 



The consequence of this is that the best fit of the data to 

m = a (K c - k) u (24) 



quoted in Jl3| (k c = 0.821(5), v = 0.89(7)) is now seen to be a poor fit to the data 
very close to the critical point. Analysis of the new data set gives 

v = 0.71(5) (25) 
k c = 0.814(2) (26) 
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with a x 2 P er degree of freedom of 1.75. This fit is shown by the middle of the three 
lines in the figure while the outer lines show the fits with v taken such that the \ 2 
per degree of freedom is 2.75. If the rh values are fitted instead then we find 

v = 0.68(10) (27) 
k c = 0.813(3) (28) 

with a x 2 P er degree of freedom of 0.98. The central values are in good agreement 
with ( p6|) but the errors are larger reflecting the larger errors in the rh values. 

It is encouraging that this new direct measurement of v is apparently in much 
better agreement with the result obtained by finite-size scaling of the specific heat 



peak 1 13, IS] but salutary in that it is necessary to work very close to the critical 



point to obtain this agreement. 



5 The specific heat 

There is a potential ambiguity in our treatment of the specific heat depending upon 
whether we calculate it with the original action S ec ([D or the action S ec (|13|) used 
in the actual simulation; we have analysed both with very similar results for critical 
exponents so here we will just present the results for S ec . The specific heat arising 
from the gaussian part of the action is analytic as a function of n so we subtract it 
and define the specific heat to be 

C = (S 2 y J - (S sym ) 2 - /3 4 ((S 2 G ) - (S G y) (29) 

Our data is shown in fig. 7. The heights and locations of the peaks for various system 



sizes can be analyzed using the standard finite size scaling relations pT 



C max = a + bL" (30) 

Kmax = K-C + b'L~~ (31) 

where a, b and b' are constants and it is expected that uj = ajv. Unfortunately our 
systems have two different shapes (this is forced upon us by our FFT implememta- 
tion which only allowed N to be a power of 2) and the coefficients b and b' might 
depend upon the shape; ignoring this for the moment and taking L = we find 



u = 0.74(20). (32) 

This can be compared with the value uj = 0.76 deduced from the symmetric lattice 
data alone so we conclude that at this level of accuracy it is justifiable to include all 
the data in order to obtain an error estimate. Although the central value for u differs 



from []T9| who found uo = 1.11(11) it is clear that with present statistics the error on 
uj is rather large. The errors on the measured values of K max are sufficiently large for 
the smaller lattice sizes that (|3TD is not very effective for deducing the correlation 
length exponent; however the data is consistent with the value v = 0.71(5) derived 
from the correlation functions. 
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Assuming the hyperscaling result a — 2(1 — u) and the finite size scaling result 
uj = ajv we find that 



v = 0.73(6) (33) 
a = 0.54(10) (34) 



Using the measured value of u and our independent measurement of v (^6]) we obtain 

a = uju = 0.53(14). (35) 

Alternatively a can be found from (|26|) by assuming hyperscaling in which case 

a = 2(1 -v) = 0.58(10). (36) 

The agreement between these methods of determination is very satisfactory but the 
value of a is much larger than that estimated by a direct fit to the specific heat 
in ||13|| . However the shape of the specific heat curve does not determine a very 
accurately and a fit of the form 

C = a + b(K c - n)- a + ck (37) 



to the full 64 2 data for « < 0.795 with k c = 0.814 fl26|) and a = 0.58 © yields a 
X 2 per degree of freedom of 1.7; the line through the 64 2 data points in fig. 7 shows 
this fit. Thus we can now conclude, on the basis of much greater statistics than in 
lfL3|| , that all the different methods of analyzing the critical behaviour are broadly 
consistent. 



6 Normal- Normal Correlation Functions 

The normal-normal correlation function defined by 

G Q (0 = (n(0 • n(0)> (38) 

where 

n(D=ti(e)xt 2 (e) (39) 

is interesting because it tells us something about the four-point couplings in the 
theory. If the model were described exactly by the free field Lagrangian d|) then 
this correlation function would be given by 



G n (0 = I (G n (0G 22 (0 - G 12 (0G 21 (0) (40) 



3 

At large distances this gives 

Corrections to this behaviour can take two generic forms. Firstly the power of the 
scaling law at large |£| may change or the coefficient of the |£| -4 behaviour could be 
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different from that predicted by (|40|); this is the behaviour that would be obtained 
if a relevant operator other than those in £f ree were present. On the other hand 
an irrelevant operator would generate deviations from (f4*0|) which fall off faster than 
|£| ~ 4 at large distances. For example, adding interaction terms so that the effective 
Lagrangian is 

Ceff = Cfree + jf^l^ * d 2 X| 2 , (42) 

where g is a dimensionless constant, gives a super-renormalizable theory and at large 
distances the induced correction to G n (£) falls off as |£| -6 . In fact our data is fully 
consistent with the picture that in the crumpled phase there are only interaction 
terms of a super-renormalizable nature. 

In fig. 8 we show the measured G n (£) for £ lying along the x-axis at k — 0.77 on 
a 64 2 lattice; the "free" line shows the form predicted by £/ ree using the parameters 
A, m determined by the fits to G\\ described in section 4. Clearly at distances 
|£| > 18 the predicted form is consistent with the measured one while at shorter 
distances there is a large discrepancy. This discrepancy can largely be explained by 
assuming a four-point coupling of the form (|42"D and computing the diagram shown 
in fig. 9. The "interacting" line in fig. 8 shows the computed G n (0 obtained in this 
way with g tuned to 4.7 which gives the best fit. As can be seen this fits the data 
much better at short distances but makes little difference to the long distance tail. 
As we move closer to the critical point it is necessary to look at larger distances 
to see the asymptotic behaviour. The behaviour of G n (£) at k — 0.79 is shown in 
fig. 10. On a 64 2 lattice there is now a discrepancy between the data and the free 
field prediction out to distances |£| ~ 25. the picture also shows what happens on 
a 128 2 lattice. Unfortunately the data is not very good but there is no evidence for 
a systematic deviation from (fiTf ) at large distances. The two curves showing the 
predictions of (f42|) on 64 2 and 128 2 lattices also show that there are significant finite 
size effects on the 64 2 lattice at this k; clearly using the 128 2 data to estimate g 
would yield a smaller value. 

The values of g deduced at some k values are shown in Table 1. These values 
should be treated with some care because it is clear that there is a systematic 
discrepancy between the data and the prediction even when the interaction term is 
included; this makes it difficult to estimate an error for the values of g quoted but 
it should be taken to be at least of order 20%. It can be seen from the table that 
the value of g is apparently falling as k approaches k c so there is no indication for 
the appearance of a relevant perturbation away from £/ ree as the system approaches 
the critical point. 

Finally we note that because the corrections to £/ ree are apparently super- 
renormalizable their presence does not affect our analysis of the two point function. 
At 0(g) only the m parameter is renormalized and the form of the two point func- 
tion is unaffected; 0(g 2 ) contributions will affect at short distances but it is 
clear that there is a substantial range of £ for which they can safely be ignored. 
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7 Discussion 



We have determined the critical exponents a and v by several different methods 
involving measurements of correlation functions and of the specific heat and making 
various theoretical assumptions. All the results obtained are in good agreement 
with each other. However in order to obtain this agreement it is necessary to make 
measurements with small statistical errors very close to the critical point and the 
nature of the analysis is slightly delicate. To obtain manifestly assumption-free 
measurements of the exponents would clearly require the use of much bigger lattices 
and need very much larger amounts of computer time than we have used. 

For the first time, and as a direct consequence of our high statistics, we have been 
able to make a convincing measurement of the normal-normal correlation function. 
The consistency of our results with the idea that the effective field theory at the 
critical point is Cf ree means that there is still an unresolved puzzle concerning these 
models. At the critical point the free energy density / on a strip of infinite length 
and finite width L with periodic boundary conditions is expected to behave as 

7T A 

If the critical point is conformally invariant then p — 1 and A is the central charge 
23|j . Attempts have been made to measure this effect |22j [T~9|1 ; assuming p = 1 these 



workers all found that A < 1. On the other hand £/ ree (which is not conformally 
invariant) predicts that p = 1 and A — 6. Whether or not the critical point is 
conformally invariant these two results are apparently inconsistent. However these 
measurements of c are all done on strips that are rather narrow, the widest strip 
being only twenty lattice spacings wide. From our measurements of G n at k = 0.79 
we know that physics differs substantially from the asymptotic behaviour implied 
by Cfree (fig- 10). Thus, if one believes the scenario implied by C e ff then at k 
values less than 0.02 away from the critical point one would not expect to see the 
correct asymptotic behaviour of ((E|). In fact these measurements are done as close 
to the critical point as possible so it may be that they are misleading. Of course, 
it is equally possible that the C e jf scenario is incorrect; our measurements are 
consistent with it but they are not accurate enough at very large distances to rule 
out other possibilities. To make further progress it will probably be necessary to 
make renormalization group studies of this model. 
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Figure Captions 



1. The value of o ([13]) required at different k values to yield p — 1 ( [To ) as 
determined on a 16 2 lattice. 

2. The ratio p (|i~5|) determined on 32 2 and 64 2 lattices. 



3. The bin size, Tun, estimated from the specific heat (17) as a function of k. 

4. The location of the zeroes of the correlation function Gn as a function of k 
for different lattice sizes. Errors are smaller than the size of the symbols. 

5. The correlation function Gn in the region where it is negative for different k 
values on a 64 2 lattice. The dashed lines show a fit of the form of (|23| ) with 
the parameters given in Table 1. Errors are discussed in the text. 

6. The mass gap m (see Table 1) plotted as a function of k. The three lines are 
discussed in the text. 

7. The specific heat as a function of k for different lattice sizes. The line is the 
fit of fl3"T|) to the 64 2 results with the constraint that k c = 0.814. 

8. The correlation function G n at k = 0.77 on a 64 2 lattice. The "free" line is 
that predicted by (jHJ) while the "interacting" line is that predicted by ( |42| ) 
with the parameters given in Table 1. 

9. Correction to G n at O(g). The open circles represent n sources, the filled circle 
a vertex of the form ( fi2|) and the lines GijS. 

10. The correlation function G n at k — 0.79 on a 64 2 and 128 2 lattices. The "free" 
line is that predicted by ([M]) while the "interacting" lines are those predicted 
by (|42|) with the parameters given in Table 1. 
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Figure 5: 
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Figure 7: 
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Figure 9: 



Figure 10: 
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